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We investigate the motion of active semifiexible filament with shape kinematics and hydrody- 
namic interaction including. Three types of filament motion are found: Translation, snaking and 
rotation. Change of flexibihty will induce instability of shape kinematics and further result in 
asymmetry of shape kinematics respect to the motion of mass center, which are responsible to a 
continuous-like transition from translation to snaking and a first-order-like transition from snaking 
to rotation, respectively. Of particular interest, we find that long-range hydrodynamic interaction 
is not necessary for filament rotation, but can enhance remarkably the parameter region for its 
appearance. This finding may provide an evidence that the experimentally found collective rotation 
of active filaments is more likely to arise from the individual property even without the long-range 
hydrodynamic interaction. 
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Due to its important role for cell division, cell motil- 
ity, force generation [l], Q and self-assembly of active soft 
matter 0, motion of filament at subcellular level has 
gained rapidly increasing attentions Usually, the 

dynamics of filament is governed by low Reynolds num- 
ber hydrodynamics, where viscous effects dominate and 
inertial effects are negligible 9] , and involves the active 
process by motor proteins walking on it, which convert 
chemical energy into work and generate propulsions. So 
far, several theoretical models have been proposed to de- 
scribe the active filament. For example, a rod model 
is used to study enhanced ordering of interacting fila- 
ment by molecular motors 6j, and a structure consists 
of a spherical head, a curved midpiece and a beating 
tail has been constructed to reveal the hydrodynam- 
ics of sperm cells near surfaces ^lOj. Lately, a lattice- 
based cellular automaton model is implemented to re- 
produce wave- like collective motion patterns [1]. Never- 
theless, these models either neglect the shape kinematics 
of each filament part evolving dynamically as a result 
of interaction between the actuation, material elasticity, 
and viscous dragQ, or miss the hydrodynamics induced 
long-range interaction both of which are important in 
understanding the active motion of filament. Recent ex- 
periments in planar geometry motility assay reveals that 
active filaments will self-organize to form coherently mov- 
ing structures, such as propagating density waves and 
vortex motion[3, H|. While Schaler et al show that long- 
range hydrodynamic interactions play a crucial role in the 
pattern forming mechanisms Q, Sumino and his collab- 
orators believe that the emergence of vortex structures 
is the result of the instinct rotation motion of single fil- 
ament in combination with local interactions and is no 
need for long-range interactions [s*]. This argument then 
calls for a complete study of motion of single active fila- 
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FIG. 1: (a) Schematic diagram of the bead-rod model of the 
active filament. 



ment. 

In this Letter, we investigate the motion of active 
filament based on a simple theoretical model includ- 
ing both long-range hydrodynamic interaction (HI) and 
shape kinematics. We find that flexibility is crucial for 
the motion of filament by controlling the shape kinemat- 
ics. For complete rigidity, activation results in a trans- 
lation movement. By a qualitative analysis and simu- 
lations, we show that, as the rigidity k decreases, fila- 
ment undergoes a transition from translation to snaking 
at K = Ki due to a flexibility induced instability of shape 
kinematics. By further decreasing rigidity, an interesting 
rotation of filament can be observed for k < K2 as a, result 
of asymmetry of shape kinematics respect to the motion 
of mass center. The scaling behavior of the continuous- 
like transition near ki is revealed, and mixed-type move- 
ment of filament is found near K2 indicating a first-order- 
like snaking-rotation transition. Particular attention is 
paid on the effect of long-range HI on the rotation of fil- 
ament. The comparison between filament motion with 
and without long-range HI shows that, long-range HI is 
not a necessary component for the filament rotation, but 
can amplify it boundary K2 considerably, which indicates 
that more rigid filament can also rotate when long-range 
HI is included. 
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Model and method-Owe simulation is based on a bead- 
rod model consists of N beads whose position is r^,? = 
1, A^, and N—1 inextensible rods with rod vector v.; — 
Ti+i — Yi,i — 1, ...^N — 1, whose schematic diagram is 
shown in Fig[T] The motion of filament is given by the 
following equation 
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where R = (ri, r^r)^. The flexibility of filament is 
realized by a bending potential 

[/■^(r„r,+i,r,+2) - = l,...,iV-2, (2) 
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where 6i is the angle between two adjacent rod vectors 
Vi-|_i and Vi. The exclusive volume effect is measured by 



0, 



rij < 2^a 
Vij > 2^/2a 



(3) 



in which rij is the distance between bead i and j. The 
activation on each rod segment is considered to be a con- 
stant force 2a along the rod vectors, which is coarsen 
from a detailed rod modelQ and neglect the transverse 
stretch of motors for simplicity. Then, the active force 
on bead i is the resultant of ones on rod i ~ 1 and i. 
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where Vi = Vi/|vi| is the unit vector of v^. HI is taken 
into account by using the Rotne-Prager-Yamakawa ap- 
proximation. The elements of mobility matrix M is given 

by 
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FIG. 2: Flexibility induced motion transition of active fila- 
ment, (a) Order paremeter Q (time-average velocity of mass 
center) as a function of flexibility k. Insert panel of (a): Three 
types of filament motion (presented by trajectory of mass cen- 
ter, solid line), translation (k — 8), snaking (k — 4) and 
rotation (k = 1.5). (b) Scaling behavior near the translation- 
snaking transition point ki . The simulation result can be well 
fitted by a line with slope 0.479 (c) A mixed-type of active 
filament motion observed near the snaking-rotation transition 
point K2. 
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The diagonal elements present the local friction and oth- 
ers the long-range interaction. The self mobility fi — 
l/(67r?7a) of a spherical particle is related to the Stokes 
friction coefficient, where a is the bead radius, and rj is 
the viscosity of the fluid. <;(t) is a set of independent (0, 1) 
Gaussian white noise, and the amplitude matrix obeys 
BB-'" — 2fcBTM. In simulation, rod constraints are usu- 
ally ensured by SHAKE algorithm^ll]. For convenience, 
we describe rods in our simulation by a FENE-fraenkel 
spring with very large spring constant [l^. 



i/(|v|-l) 



l-(l-|v|)2/a2 |v| 



,(l-a)<|v|<(l+a), (6) 



where H is the spring constant and a is the extensibility 
parameter that defines the maximum possible deviation 
between the actual spring length and the natural length. 



To prevent overstretching or overcompressing the springs, 
an implicit predictor-corrector method is employed to 
calculate the rod length, which has been described in de- 
tail in[l^. The length, time and energy scale are chosen 
as diameter of bead 2a, integration time step dt and ^. 
We fix iV = 10, L = 2{N - 1) = 18 and a = 4 if not oth- 
erwise stated. All the results are derived by simulations 
for at least 10'' time steps. 

Result-Three types of motions (presented by trajec- 
tory of mass center of the filament) are observed in differ- 
ent region of filament flexibility. In the right insert panel 
of Fig[2l^a), the filament is of straight line shape, and also 
moves forward along a straight line (translation). The 
middle insert panel of Figl^Ja) shows a snaking move- 
ment of the filament while the shape of filament is no 
longer a straight line but becomes cilialike beating. In- 
teresting motion is found when the filament is flexible 
enough (the left insert panel of Figl^fa)). The filament 
rotates spontaneously within the effect of activation, flex- 
ibility and HI. 

To take a global picture of the filament motion, we 
define an order parameter to distinct these three types of 
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movements. 

Q = — I lim / RMcdt\, (7) 

where Rmc is the coordinate of mass center, and Cq is 
value of stationary rate of translation movement for the 
same active force. In Figl^Ja), we present a quantita- 
tive description of the filament motion as a function of 
flexibility. It is found that, as k decreases, filament un- 
dergoes two motion transitions: the first transition from 
translation to snaking at ki = 6.0, and the second from 
snaking to rotation at K2 — 1.9. A scaling behavior 
near the continuous-like transition point ki is given in 
FiglUb). The simulated dependence of log\{Q — Qi)/Qi\ 
on log\{K — ki)/ki\ can be well fitted by a line with slope 
0.479, implying that the snaking movement may arises 
from an instability of translation occuring via a super- 
critical bifurcation. In Fig|2Jc), we record a mixed- type 
of motion near K2. The filament first snakes for a while, 
then rotates for a short period and turns back to snaking. 
This process is repeated over and over again. The coexis- 
tence of snaking and rotation indicates that the filament 
motion is multistable near the transition point K2 , which 
also suggests that the transition from snaking to rotation 
may be a first-order like transition. 

We then do further studies about the underlying mech- 
anism of these two transitions. The result reveals that 
motion transitions are closely related with the shape 
kinematics of filament. Firstly, a qualitative linear sta- 
bility analysis can be given by taking the last 3-bead 
segment as a sample. When the shape of straight line is 
stable and neglect long-range HI, the shape kinematics 
of the 3-bead segment is given by oc asinO — kO/tt. The 
straight line solution 9 = loses its stability for k < an. 
To verify this instability, the distributions of the angle ip 
between velocity of the first bead and of the mass center 
are shown in Figl^l^a). For translation, if only distributes 
near ip = 0. Once filament moves in type of snaking, 
two symmetric peaks of ip emerge, which indicates that 
straight line shape does lose its stability at k = ki, and 
results in a shape kinematics of cilialike beating. The 
distributions for k = 3, 4 and 5 show that decreasing of 
K amplifies the amplitude of cilialike beating. Then, we 
then simulate the distribution of p for rotation (FiglJl^c)). 
Similar to the one of snaking, the distribution shows two 
peaks, too, and the amplitude of cilialike beating is fur- 
ther enlarged. However, these two peaks are no longer 
symmetric to each other in both positions and the rela- 
tive probabilities. As a result of the asymmetry of shape 
kinematics respect to the motion of mass center, the fil- 
ament keeps tuning clockwise or anticlockwise which is 
determined by the preference of the distribution. Cumu- 
lation of the tuning finally forms the rotation motion of 
active filament. 

Now, we try to reveal the effect of long-range HI on 
the motion of active filament. In our simulation, long- 
range HI can be conveniently turned off' by considering 
only the diagonal elements of mobility matrix M. Sur- 
prisingly, in contrary to the intuition that long-range HI 




-1.5 -1.0 -0.6 0.0 0.5 1.0 

9 



FIG. 3: The stability and symmetry of shape kinematics re- 
spect to the motion of mass center, which is presented by the 
distribution of angle p between velocity of the first bead and 
the mass center for the motion of (a) snaking and (b) rota- 
tion. The distribution of p for translation is also plotted in 
(a) in patterned bar. 



plays important role in formation of different types of fil- 
ament motion, all the three types of filament motion are 
still observed in simulation without long-range HI. The 
transitions from translation to snaking and from snaking 
to rotation are similar to their respective one with long- 
range HI, too. To show this result clearly, we plot the 
phase diagram of active filament motion with and with- 
out long-range HI in a — k plane in FigUJa). It can be 
seen that, no matter long-range HI presents or not, three 
distinct motion regions can be identified clearly. In other 
word, long-range HI is not relevant to the formation of 
active filament motion. Besides of this, several other ob- 
servation can be concluded from Fig|4lja). Firstly, the 
instability of straight line shape occurs at ki oc a, which 
is accordant with the qualitative linear stability analy- 
sis. Determining of the exact instability condition may 
be mathematically complicated, and needs further stud- 
ies. Secondly, both the transition boundaries with and 
without long-range HI are proportional to the active force 
except that the coefficient are different. This implies that 
the long-range HI will result in an effective coefficient of 
the instability. Finally, the filament motion in type of 
rotation or snaking are enhanced by long-range HI. By a 
comparison of flexibility region in which filament moves 
in type of rotation or snaking with and without long- 
range HI for different active force a, it can be seen that, 
when long-range HI presents, the range of rotation re- 
gion K < K2 is of magnitude larger than the one only 
considering local HI. In the meanwhile, width of snaking 
region ki — K2 is also enlarged several times by long-range 
HI. Thus, FiglU^a) gives evidences that single filament 
does rotate without long-range HI, and the parameter 
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FIG. 4: Phase diagram of active filament motion on (a) a — k 
and (b) L — k plane without (left) and with (right) long-range 
HI (LRHI). T: Translation. S: Snaking. R: Rotation. The 
translation-snaking (T-S) and snaking-rotation (S-R) bound- 
ary are plotted by line-cycle and line-square respectively. We 
fix filament length L = 18 for (a) and active force a = 4 for 
(b). 

region for rotation can be enhanced by long-range HI. 
This finding supports that the shape kinematics is suffi- 
cient to form instinct rotation of single active filament, 
and long-range HI may be not necessary for the forma- 
tion of collective vortex of filamentsfsl]. 

Besides of long-range HI, the effect of activation a 
on the motion of filament can also be concluded from 
FigHfa). For a given k, a will induce motion transition of 
filament, too. Similar to decreasing of k, the filament first 
undergoes transition from translation to snaking, then 
the transition from snaking to rotation, by increasing a. 
It is noted that, when long-range HI absents, the slope of 
snaking-rotation boundary is very small (FigUJ^a)). As a 
result, rotation may not be observed by finite increasing 
of a. When long-range HI presents, however, the slope of 
snaking-rotation boundary becomes much sharper, and 



rotation will be more easily observed in simulation. At 
last, the effect of filament length L is studied. The sim- 
ulated phase diagram in L — k plane in FigUfb) shows 
that, when long-range HI presents, proportionable ki is 
needed for the transition from translation to snaking as a 
result of increasing of L. The snaking-rotation transition 
point K2, in contrary, saturates at about K2 = 2.6. Com- 
pared with ones considering only local HI, long-range HI 
can also enhances the rotation and snaking region. In 
addition, we also simulate the motion of much longer fil- 
ament. When K is small, the filament rolls into a coil 
state, and rotation is hardly to be observed in simula- 
tions no matter including long-range HI or not. On the 
other hand, for sufficient large k, the translation is still 
stable and transition from translation to snaking can al- 
ways present. 

Conclusion and discMssion-Flexibility induced motion 
transition has been studied by a simple active filament 
model which allows consideration of both long-range HI 
and shape kinematics. This research reveals that the 
motion of filament is mainly determined by flexibility via 
changing the stability and symmetry of shape kinemat- 
ics. We also demonstrate that long-range HI can en- 
hance but not crucial for the motion of single filament 
such as rotation. These findings may be helpful to un- 
cover the underlying mechanism of collective behaviors 
observed in experiments. What's more, based on the sim- 
ple model of single active filament, building up a step by 
step theoretical framework which relates the individual 
shape kinematics and long-range HI to the nonequilib- 
rium self-assembly of active filaments is possible. 

When we finish this Letter, we notice that similar work 
on active filament has been done by Jayaraman et al and 
reported in the latest Physical Review Letters|lJ|. They 
have found that the interplay between long-range hydro- 
dynamics and semiflexibility is necessary for the rotation 
motion of the filament: A bending instability of an ini- 
tially straight filament spontaneously breaks flow sym- 
metry and leads to autonomous filament motion which, 
depending on conformational symmetry, can be transla- 
tional or rotational. The difference role of long-range HI 
revealed in their work and ours may need further under- 
standing. 
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